Consequences of extended s±-wave pairing in iron-based superconductors 
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Motivated by a recent experiment of Song et al. [Science 332, 1410 (2011)], we theoretically study 
the spin dynamics, charge dynamics, and point-contact Andreev-reflection spectroscopy (PCARS) 
of two-band iron-based superconductors of a possible extended s±-wave pairing symmetry. We 
consider the case of a dominant s± gap blended by a secondary extended s component in which gap 
nodes can develop in the Fermi pockets near zone corner and/or boundary. Due to the strong nesting 
effect associated with nodal regions, dynamical spin and charge susceptibilities can exhibit strong 
peaks at momenta near (±7r/2,0), (±7r,±7r/2), as well as (±7r,0) in the unfolded Brillouin zone. 
For PCARS, considering an anisotropic band effect induced by an applied voltage, [100] differential 
conductance can exhibit a V-shape behavior manifesting a gap node occurring in such direction. It 
is highly suggested that the above features can be experimentally investigated to help sorting out 
the pairing symmetry of iron-based superconductors. 

PACS numbers: 74.25. Ha, 74.20. Mn, 74.20.Rp, 74.50. +r 
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I. INTRODUCTION 

Pairing symmetry of the Fe-based pnictide and chalco- 
genide superconductors is currently a hot topic in the 
study of superconductivity. While more and more ex- 
perimental results have suggested that the order param- 
eter in these materials is likely to be fully gapped s±- 
wave, 1-4 whether there is a node in the gap remains 
controversially 5,6 For instance, in a recent scanning tun- 
neling microscopy (STM) measurement of Song et al., 7 a 
nodal and two-fold symmetry gap is revealed in an iron 
selenide (FeSe) superconductor. It is no doubt that a 
central issue towards understanding the iron-based su- 
perconductors (FeSCs) is to unambiguously identify the 
pairing symmetry of these materials. 

Among many different probes, inelastic neutron scat- 
tering (INS) measures the two-particle excitations and 
can give direct information on the momentum and en- 
ergy dependence of the quasiparticle excitation and the 
pairing gap. A strong coherence peak can emerge in 
the dynamic spin susceptibility if the corresponding two- 
particle excitation is highly degenerate (i.e., in good nest- 
ing condition). In addition, strong peaks in INS can oc- 
cur due to the resonant nature. A conclusive INS mea- 
surement of FeSCs remains unsettled however to which 
wave vectors of the strong peaks are observed to be re- 
markably material dependent. Some INS measurements 
have reported that spin resonances occur at the wave 
vector near (it, tt) [or (tt, 0)] in the folded (or unfolded) 
Brillouin zone (BZ) in llll, 8 122, 9 ' 10 111, 11 and 11 
families. 12 ' 13 Other INS measurements on the new 122* 
family Aa;Fe2Se2 (A=K, Rb, and Cs) 14 of electron Fermi 
surfaces only 15 ' 16 have revealed that a resonance peak 
occurs at the wave vector (tt, tt/2) in the unfolded BZ. 17 
On the theoretical side, on the other hand, prediction of 
the positions of the resonance peak in momentum space 
remains controversial. It has been predicted that for s± 
pairing state, a strong coherence peak can exist in the 



dynamic spin susceptibility x(Q, u ) a t the nesting wave 
vector Q=(0,tt) or (tt,0) in the unfolded BZ. 18-20 How- 
ever, based on a pairing potential associated with the 
predicted x(Q, u>), the SC state are found to be incon- 
sistent with the previous prediction. 21,22 

Charge dynamics which is accessible by high-energy 
electron scattering or X-ray scattering is another ideal 
candidate for studying the pairing symmetry of iron- 
based superconductors. While charge and spin suscep- 
tibilities are coupled to different coherence factors, as far 
as two-particle excitation is concerned they are qualita- 
tively similar. Thus for a comparison point of view, it 
is also useful to theoretically study the charge dynamics 
within the same framework. 

Another high-resolution phase-sensitive probe to de- 
tect the pairing symmetry is the point-contact Andreev- 
rcflcction spectroscopy (PCARS). 23-25 The situations is 
still too early to make a conclusion however. Some 
PCARS measurements showed two coherent peaks and 
indicated that SC pairing state might be fully gaped on 
the Fermi surface (FS). 26-28 Others showed a zero-bias 
conductance peak (ZBCP) and implied the presence of 
zero-energy bound states or Andreev bound state (ABS) 
on the interface. 29-31 Moreover, depending on the direc- 
tion of the sample interface, some PCARSs have shown 
ZBCP coexisting with finite-energy coherent peaks. 32,33 
Most of theoretical studies so far have focused on the ex- 
planation of the ZBCP. 34,35 It is important to carefully 
identify whether a gap node exists through the PCARS 
data. 

In the current paper, we use a minimal two-orbital 
model 36 to study the spin dynamics, charge dynamics, 
and PCARS of FeSCs. Motivated by the recent STM 
experiment 7 mentioned previously, we consider a pair- 
ing gap of a primary s±-wave component plus a sec- 
ondary extended s-wave component, called the extended 
s_j_-wave state. In fact, this mixed pairing state is sup- 
ported by a theoretical work of Yang et al. 37 who did a 
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variational quantum Monte Carlo calculation and con- 
cluded that both s-|--wave and extended s-wave pair- 
ings are equally energetically favorable in the five-band 
FeSCs. In the studies of spin and charge dynamics, it 
will be shown that in addition to wave vectors (±7r, 0) 
that most previous works focused on, strong coherence 
peaks can also occur around wave vectors (±7r/2,0) and 
(±7r, ±7r/2), which is a unique feature in the extended 
s±-wave state. In the studies of PCARS, considering an 
external anisotropy effect due to an applied bias voltage 
that leads to a relatively larger (smaller) Fermi pocket 
for Pi- (/?2-) band and hence a node can develop in the 
(it, 0) direction in /?i-band, 38 low-energy differential con- 
ductance dl/dV along the [100] direction will feature a 
V~-shape curve and manifests the existence of a node. 

This paper is organized as follows. In Sec. II, we in- 
troduce the two-orbital model and especially show how 
the mixed gap behaves in different FSs as the secondary 
extended s-wave component changes. Sec. Ill is devoted 
to study the dynamical spin and charge susceptibilities 
for the model introduced in Sec. II. Section IV gives a 
theoretical study of the PCARS. Sec. V is a brief sum- 
mary. For self sustainability, Appendix A gives detailed 
forms of the irreducible spin and charge response func- 
tions of a two-band superconductor. A brief discussion 
of Random-Phase Approximation (RPA) on the vertex- 
corrected spin and charge response functions is given in 
Appendix B. 



II. MODEL 

We consider a minimal two-orbital model for iron- 
based superconductors in which both d xz and d yz or- 
bitals, coupled by the d xy orbital, arc considered in a 
two-dimensional square lattice. The Hamiltonian is 36 



e x (k) - fi eay(k) 
e xy (k) e„(k) - /i 
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where ^ [dL(k), 4 f (k) ] with 4 ff (k) [d^k)} ere- 
ating an electron in orbital d xz (d yz ) of wave vector k 
and spin a and 
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Here t\, £2 are the nearest-neighbor hoppings and is, £4 
are the next-nearest-neighbor hoppings. After Bogoli- 
ubov transformation, Hamiltonian (1) becomes 



H Q = J2 e,(k) 7 L(k)7-(k), 



k.(T,i/— ± 

where the band dispersions 
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FIG. 1. (Color online) Frame (a): Fermi surfaces of the two- 
band model in the unfolded BZ. Red (blue) curves correspond 
to the hole (electron) Fermi pockets Qi and a.2 (Pi and /fe). 
Frame (b)-(d): Gap magnitudes around the 02 (red) and Pi 
(blue) FSs for mixing ratio (b) x = 0, (c) x = 0.2, and (d) 
x = 0.4 respectively. 

with e±(k) = [e x (k) ± e y (k)}/2 and jua(k) is the band 
operator which destroys a quasiparticle of wave vector 
k and spin a in band v. Operators d ra (k) (r = x,y) 
and 7^0- (k) (v = ±) are related to each other via the 
transformation 



d ra (k) = <(k)7^(k), 



(5) 



u=± 



where 



^k) = a^k)=sgn[ e;c!/ (k)] 



.(k) 
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Considering a spin-singlet pairing state, the following 
BCS mean-filed interaction Hamiltonian 



H SC = - J2 [A,(k) 7 t t (k) 7 ^(-k)+H, 



(7) 



k,i/=± 



can be added to Hq in (3). We consider the gap function 
being in the extended s± state: 

A, 

A^(k) = Ai cos k x cos k y + — (cos k x + cos k y ) 



An 



(1 — x) cos cos ky + — (cosfc x + cos k y ) (8) 
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where a secondary extended s component (oc A2) is 
added to the dominant s ± component (oc A!). 7,37 A 
"mixing ratio" x is introduced in the second line of (8) for 
convenience. For simplicity, the gap functions are taken 
to be the same in both bands. 

Fig. 1(a) plots the FSs of the two-band model in which 
there are two hole Fermi pockets [a± and «2 associated 
with £_(fcp) = 0] and two electron Fermi pockets and 
02 associated with £ + (fc F ) = 0]. The parameters used are 
t\ = —0.8, ti = 1.4, £3 = < 4 = —0.9, and chemical po- 
tential [i = 1.15. 38 All energies are in units of t(> 0) 
which is material dependent. In Figs. l(b)-l(d), we plot 
the gap magnitudes associated with 012- and /3i-band FSs 
by varying the mixing ratio x and with a fixed ampli- 
tude A = 0.05 (also in units of t). The gaps associated 
with oil- and /3i-band FSs are of particular interest in 
the current context, which will become clear shortly. As 
shown in Figs. l(b)-l(d), gap magnitude on the a-i FS 
decreases as the mixing ratio x increases. In particular, 
in the case of x — 0.4 [see Fig. 1(d)], the gap magnitude 
around a 2 band can reduce to zero {i.e., becomes nodal). 
Note that the gap magnitude around the ot\ band (not 
shown) remains roughly unchanged as x changes. In ad- 
dition, of equal importance, the gap magnitude around 
the /3i-band FS becomes smaller and smaller in the k x 
axis as x increases. When x is increased to be roughly 
0.7, or when there is an induced band anisotropy due to 
an applied voltage with x = 0.4 (see Fig. 5), a gap node 
can actually develop in the /?i Fermi pocket in the k x 
direction. 



III. SPIN AND CHARGE DYNAMICS 

This section is devoted to study the dynamical spin and 
charge susceptibilities of FeSCs. To save the space here 
and still for self sustainability, detailed forms of the irre- 
ducible spin and charge response functions of a two-band 
superconductor are to be given in Appendix A. A brief 
discussion of the Random-Phase Approximation (RPA) 
on the vertex-corrected spin and charge response func- 
tions is given in Appendix B. 



A. Dynamic Spin Susceptibility 

We first consider the dynamical spin susceptibility in 
the T — > limit. For spin-singlet pairing, dynamical spin 
susceptibility is proportional to the imaginary part of the 
spin response function defined as 

i-P 

X + _(q,ia;„) = V / dre^ r (r T 5+(q,r)^-(-q,0)), 
r,t Jo 

(9) 
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FIG. 2. (Color online) (a) Irreducible dynamic spin sus- 
ceptibility ImXs(q,<^ = 0.3Ao) for wavevectors q along the 
(0, 0) ->■ (tt, 0) (71- , vr) ->■ (0, 0) route in the BZ and for 
different mixing ratio x. (b) Contour plot of the excitation 
spectrum — oj/2 = 0.15Ao with x = 0.4. Blue arrows 
(P2-P4) correspond to nesting wave vectors between the hole 
o?2 (black) and electron /3i (red) bands. The one green arrow 
(pi) corresponds to a nesting wave vector within the electron 
fix band. These vectors correspond to those peaks emerging 
in frame (a). 

where r,t = x,y are the orbital indices and the spin op- 
erator 

S+(q, t) = J2 4 t ( k + <*> T Ki( k > T ) 
k 

= [5-(-q,r)]t. (10) 

For simplicity, we denote Xs = Xh throughout this pa- 
per. In the T — > limit, only the first and second terms 
contribute in the irreducible spin response function x° 
given in Eq. (Al) in Appendix A. Furthermore, only the 
first term contributes if uj > is concerned. Conse- 
quently, apart from the effect of the coherence factors, 
the contributions to Imx° are dominated by the condi- 
tion of the particle-particle excitation 

cj = ^(k + q)+E„(k). (11) 

Here E v (k) = ^^{k) + |A„(k)| 2 is the quasiparticle ex- 
citation spectrum with the band indix v = +, — . 
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FIG. 3. (Color online) Comparison of the spectra in irre- FIG. 4. (Color online) Comparison of the spectra in irre- 
ducible and RPA-corrected spin susceptibility ImXs and Im^s ducible spin and charge susceptibilities Im^s an d Im^c for 
for x = 0.4 and u) = 0.3Ao. Coulomb repulsion is taken to be x = 0.4 as cj = 0.3Aq. 
17 = 1.8 for Imxs- 



In Fig. 2(a), we show the calculating results of the ir- 
reducible spin susceptibility Im^s for four different mix- 
ing ratio x and a fixed (low) energy u = 0.3Ao. This 
two-dimensional plot is made by connecting important 
symmetry points in the BZ. As shown in Fig. 2(a), Imx^ 
has strong peaks emerging at wave vectors near (jr/2,0), 
(tt, 0), and (tt,tt/2). Of most interest, with the same 
broadening S = 0.01 (in units of t) and the same gap am- 
plitude A used in all cases, the peaks are seen to be much 
stronger in the larger x cases compared to those in the 
pure s±-wave (x = 0) case. To understand these phenom- 
ena, one just considers the low-energy particle-particle 
excitations with respect to the gap mixing. For pure 
s±-wave pairing (x = 0) case, Fermi surfaces are fully 
gapped and consequently low-energy particle-particle ex- 
citation is suppressed. In contrast, for extended s±-wave 
pairing of larger x, gap nodes can develop around the 
«2-band FSs as well as in certain k areas around the 
/3-band FSs. Consequently low-energy particle-particle 
excitation is enabled. (Note that no low-energy excita- 
tion with lo < Ao is allowed in connection with the fully 
gapped ai-band FS regardless of what the mixing ratio 
x is.) 

To understand the origins of the peaks further, in 
Fig. 2(b) we plot energy contours of the quasiparticle 
excitation = oj/2 = 0.15Ao with x = 0.4. The ar- 
rows denote the nesting wave vectors for which degen- 
erate strong two-particle excitations are associated with. 
These particular wave vectors, labeled by p\ ~ p±, are ex- 
actly those that the strong peaks are associated with in 
Fig. 2(a). As a matter of the fact, strong peaks emerge 
due to the scattering between electron and hole Fermi 
pockets (blue vectors) as well as within the same elec- 
tron Fermi pocket (green vector) . 

To see how the behaviors of the full spin susceptibil- 
ity differ from those of the irreducible one, in Fig. 3 we 
compare the results of Inix!? and Im^s for x = 0.4 and 
lo = 0.3A . Here the full spin response function \ s is cal- 
culated under RPA. A brief discussion on how Xs is ob- 



tained is given in Appendix B. In calculating Xsi we have 
set the intro-orbital Coulomb repulsion U = 1.8 (in units 
of t) and the Hund's coupling J = (see Appendix B). In 
view of Fig. 3, while enhancement occurs for Inr\; s , the 
results are qualitatively similar between Im^s and Im^s- 
Concerning the results of Im^s and/or Imx s in Fig. 3, 
we note the following two facts. First, the spectra for a 
higher (lower) energy with u> > 0.3A (to < 0.3A ) are 
basically the same as those of the w = 0.3Ao case. The 
only major difference is that the intensity of the peaks 
is stronger (weaker) in the higher (lower) energy case. 
Secondly, in the parameter regime we are studying, the 
peak around q = (jr, 0) is enhanced most in Imxs- Per- 
haps this is the case of most interest as most current 
INS experiments have focused on the possible magnetic 
peaks around q = (tt, 0). 8 ~ 13 We emphasize that strong 
magnetic peaks can also emerge at wave vectors around 
(tt/2, 0) and (tt, tt/2) if an extended s±-wave pairing state 
of a larger mixing ratio x is present. 

B. Dynamic Charge Susceptibility 

We now turn to study the dynamical charge suscepti- 
bility in the T — > limit. Dynamical charge susceptibil- 
ity is proportional to the imaginary part of the charge 
response function 

Xc(q,;0 = V / dTC^ r (T T p r (q,T)p t (-q,0)),(12) 

r,t J 

where the density operator 

Pr(q,r) = ^4 <T (k + q,T)d r(T (k,T). (13) 

kcr 

The detailed form of the irreducible charge response func- 
tion, x°, is given in Eq. (A2) in Appendix A. In Fig. 4, we 
show and compare the spectra of the irreducible charge 
and spin susceptibilities Im^c and Imx^ for x = 0.4 and 
u> = 0.3Aq. One sees clearly that similar strong peaks 
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emerge at wave vectors near (w/2, 0), (-7T, 7r/2) as well 
as (n, 0) in both Imx° and Imx^. Therefore the charge 
peaks with wave vectors near (ir/2, 0), (tt, 7r/2), and (7r, 0) 
can also be considered as the signature for the extended 
s±-wave pairing state. One is noted a key difference in 
the spectra of Imx° and Imx° however. In ImXs, the 
strongest peak is associated with wavevector pi , while in 
Im^c, the strongest peak is associated with wavevector 
Pi. The difference is due to different coherence factors 
coupled to Xc an d X° [ see t ne nrs t term of x° and Xc m 
Eqs. (Al) and (A2)]. 



IV. POINT-CONTACT 
ANDREEV- REFLECTION SPECTROSCOPY 




FIG. 5. (Color online) (a) Fermi surfaces of the band disper- 
sions given in (17). The parameters are the same as those in 
Fig. 1(a) together with an anisotropy Se — 0.3. (b) Gap mag- 
nitudes of the extended s±-wave gap (8) associated with Q2- 
(red) and /Vband (blue) FSs with x = 0.4 and A = 0.05. 



When a free electron is injected from a normal metal 
into a superconductor (forming a interface at x = 0) with 
an incident angle 8, there are two possible reflections. 
One is the normal reflection (reflected as electrons) and 
another is the Andreev reflection (reflected as holes, due 
to electron and hole coupling in the k subspace). The 
resulting wave function vf^ in the normal side (x < 0) 
can be obtained by the superposition of the above two 
kinds of reflected waves and the incident wave. In the 
superconducting side (x > 0), the resulting wave function 
^5 is the superposition of the electron-like and hole-like 
quasiparticle wave functions, which can be obtained by 
solving the Andreev equation. 

On key issue for the superconductor is that electron- 
like and hole-like quasiparticles will experience the effec- 
tive pairing potential A(0) and A(n — 0) respectively. 
As a matter of the fact, resulting tunneling conductance 
will depend strongly on the pairing symmetry and the 
direction of the tunneling current. Considering tunnel- 
ing current along the [100] direction, A{0) = A(ir — 0) 
for the current extended s±-wavc pairing and tunneling 
conductance will feature a nodeless gap if the mixing ra- 
tio x is small. However, as mentioned before, when x is 
large or when there is a large band anisotropy occurring 
due to an applied voltage (see later), the gap associated 
with the /3i-band will shrink in the k x direction and could 
eventually develop a node in such direction. This nodal 
gap will then be manifested in the tunneling conductance 
along such direction. 

By matching the normal and superconducting wave 
functions and their derivatives at the interface x = 0: 



ipN (r) | x =o-= ips (r) U =0 +, 

2mH #sW, 

-ips (r) U=o+= — — — |x=o+ 



(14) 



w 2 



dx 



dipN (r ) , 
dx 



\x=0- 



where H denotes the strength of a delta- function like bar- 
rier potential V(x) = HS(x), one can solve the normal 
and Andreev reflection coefficients and ta- With the 
solved rjv and and based on the well-known formula 
for the conductance 0-5 = 1 + \ta\ 2 — |fjv| 2 (see, for ex- 
ample, Ref. [39]), the normalized differential conductance 



dl/dV can be obtained to be 

16(1 + T 2 ) cos 4 + AZ 2 (1 - T 2 ) cos 2 



a s (E, 1 



|4cos 2 (9 + Z 2 (l -r 2 )| 2 



(15) 

where Z = 2mH/ti 2 is the effective potential barrier and 



E 



E 



|A(0)| 



|A(0)| 



(16) 



In a real experimental setup, the above derivation is 
equivalent to consider an STM tip made of a normal 
metal and being well point-contacted to the supercon- 
ductor. Similar derivations of Eqs. (15) and (16) can be 
found in Ref. [24] where a single-band system was con- 
sidered. 

In the current context, for simplicity, we assume that 
the two bands are decoupled completely. Nevertheless 
the resulting tunneling conductance is considered to be 
contributed by all ot\, a.2, Pi, and @2 FSs. To see the pos- 
sible nodal effect for the conductance current along the 
[100] direction, we consider a band anisotropy effect in- 
duced by an applied voltage. The applied voltage should 
be much larger than the pairing gap magnitude. 38 Con- 
sequently the induced anisotropic band dispersion can be 
approximated by 



£±(k) = e+ (k)±,/(e-(k)+<5e) 



6 2 ,(k)-/l, (17) 



where Se denotes the induced anisotropy. In Fig. 5(a), 
FSs of the anisotropic band (17) are shown in which Ri- 
band Fermi pocket becomes larger, while /32-band Fermi 
pocket becomes smaller as compared to those in Fig. 1(a). 
The parameters used in Fig. 5(a) are the same as those 
in Fig. 1(a), i.e., ti = -0.8, t 2 = 1.4, t 3 = U = -0.9, and 
\i = 1.15. In addition, the induced anisotropy is chosen to 
be Se = 0.3. Fig. 5(b) shows how the extended s±-wave 
gap behaves in the a.2- and /?i-band FSs with x = 0.4 
and Ao = 0.05. The most important feature in these 
gaps is that gap nodes develop in both a 2 - an d /3i-band 
FSs. For the a2-band gap, the node appears to be nodal 
almost everywhere around the ui FS. In contrast, for the 
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FIG. 6. (Color online) (a) Resulting [100] differential conduc- 
tance dl/dV vs. tunneling voltage E (black line). Parameters 
are the same as those in Fig. 5. For comparison, contribu- 
tions due to four Fermi pockets are plotted respectively in 
which /?2 contribution is magnified by a factor of 20 for a bet- 
ter view, (b) Comparison of resulting [100] dI/dV(E) curves 
with x — 0.4, 0.2, and 0. For clarity the curve of x = 0.2 (0) 
is shifted by 0.02 (0.04) vertically. 



ft-band gap, the node appears in the k x direction only 
which should have a strong effect on the conductance 
current along the [100] direction. 

Figure 6(a) plots the [100] differential conductance 
dl/dV versus the tunneling voltage E for a normal- 
superconducting junction. The superconductor is made 
of the two-band iron-based with the induced anisotropic 
band dispersions given in (17) and the extended s±-wave 
pairing gap given in (8). The parameters used are the 
same as those in Fig. 5. In addition, the strength of bar- 
rier is taken to be H = 10 (in units of t). To be more ex- 
plicitly, dI/dV(E) is obtained by integrating <jg(E,8) in 
(15) over the 9 angle from — tt/2 to n/2. For comparison 
purpose, we have plotted the contribution respectively 
due to ai, 012, ft, and ft bands. Of most interest is that, 
in the low E regime, in contrast to ai, oti-, and ft bands 
whose contribution is either little or almost independent 
of E, the contribution due to ft band is dominant and 
linear in E. This can be understood in the following, 
ft-band is located far away from the k x axis, thus its 
contribution (green line) to the differential conductance 
is the least. For the a^-band, the gap is nodal almost ev- 
erywhere on the FS, so its contribution to dI/dV(E) (red 
line) is almost independent of E (i.e., obeying the Ohm's 
law). For the ai-band, finite s-wave gap is fully devel- 
oped in its FS, consequently its corresponding dI/dV(E) 
(blue line) exhibits a U shape. In contrast for the ft- 
band, gap node develops in the k x direction, so for the 
[100] differential conductance, it exhibits a V shape (pur- 
ple line) and manifests the existence of a gap node along 
such direction. The overall U-shape dI/dV(E) at low E 
seems to be in consistence with the STM measurement 
reported in Ref. [7]. 

In Fig. 6(b), we show the [100] dl/dV versus E by 
changing the gap mixing ratio x. It is evident that when 



x = for a pure s±-wave pairing, all FSs are fully gapped 
and consequently the resulting dI/dV(E) exhibits a U- 
shapc behavior. Those tiny peaks corresponds to quasi- 
particle excitations associated with different Fermi sur- 
faces. It is in great contrast to the case of x = 0.4 in 
which dI/dV(E) exhibits a V-shape behavior, as already 
shown in Fig. 6(a). 



V. SUMMARY 

In this paper, motivated by the recent experiment 
of Song et al7 , we intend to study the spin dynamics, 
charge dynamics, and point-contact Andreev-reflection 
spectroscopy (PCARS) of iron-based superconductors of 
a possible extended s±-wave pairing symmetry. We con- 
sider a minimal two-band model in which superconduct- 
ing gap is dominated by the s± component but mixed by 
a secondary extended s component. Comparing to the 
pure s±-wave gap of finite gaps around all four Fermi 
pockets, in the present extended s±-wave pairing state, 
gap nodes can develop in the Fermi pockets near the zone 
corner and/or the zone boundary. Useful signatures as- 
sociated with the possible nodes are identified in spin 
dynamics, charge dynamics, as well as in the PCARS. It 
is hoped that the features discussed in this paper can be 
carefully investigated by experiments to help sorting out 
the pairing symmetry of iron-based superconductors. 

Finally we remark that some recent inelastic neutron 
scattering experiments on 1^.1^02802 samples have uncov- 
ered magnetic resonance peak at q = (tt, it/ 2) in unfolded 
BZ. 17 This implies that the pairing state in these materi- 
als is likely to be d-wave. Fermi surface topology in these 
materials is, in fact, quite different from the iron-based 
superconductors that we are considering. 
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Appendix A: Irreducible Spin and Charge Response 
Functions 



Within the two-orbital model proposed by Raghu et 
al., 36 irreducible spin and charge response functions can 
be derived to be (see also the derivation in Ref. [40]) 



7 



X?(q,-) = -^EE E ^T'(k,q) 

k r,t i/,i/'=± 

/(^(k))-/(^(k + q)) 



.4 



l-/(£„(k))-/(£„(k + q)) 



-Br 



l-f(E u ,(k))-f(E u (k + q)) 



- C, 



kq w - (k) + EvQt + q) + i(S 



kq w - ^(k) - £„(k + q) + id k -i w + ^(k) + E v (k + c 

ExLt(q^) (Al) 



f(E v ,(k))-f(E v (k + q)) 



£„<(k) -E v {k + q) + i8 



x°(q,w) 



^EE E ^T'(k>q) a 

_ f(E u ,(k))-f(E u (k + q)) 



l-f(E u ,(k))-f(E v (k + q)) 



k,q" 



- B 



k ^o;-^(k)+^(k + q) 



i<5 



D 



E v ,(k)-E v (k + q)+i6 
_ f(E v ,(k))-f(E v (k + q)) 



k,q" 



k «u + E l/ ,{k)-E v (k + < l ) 



iS 



l-/(^(k))-/(^(k + q)) 
■jj + E„i (k) + E v (k + q) + i<5 

Ex°,, t (q,^), (A2) 



where r, t = x,y correspond to the orbital indices, 
v,v' = ± correspond to the band indices, TV is 
the renormalization factor, and 5 is the broadening. 
f{E v (k)) is the Fermi distribution function with E u (k) = 
y/££(k) + |A(k)| 2 the quasiparticle excitation spectrum, 
Ci/(k) the band dispersion, and A„(k) the superconduct- 
ing gap. The coherence factors in (Al) and (A2) are 
respectively 



A ± = - 



± 



1 



G'(k) 
A I/ /(k)A I/ (k + q 



1 - 



k,q 



£„/(k)£7„(k + q. 



± 



A^(k)A„(k + q 



E„,(k)E u {k + q) 
, C'(k) 



k.q 



± 



S„/(k) 

A^(k)A„(k + q 



^(k)^(k + q) 

e,'(k) 



k,q 



1 - 



± 



E v ,(k) 
A^(k)A y (k + q 



1 - 



£„(k4 


-q) 


£„(k- 




6(M 


-q) 


£„(k- 


Hq) 




-q) 


E„{k- 


Hq) 




-q) 


£*(k- 


Fq) 



(A3) 



£^(k)£„(k + q 
and the function 

(k, q) = <(k + q)<, (k)ot, (k)ot(k + q), (A4) 
with a£(k) defined in Eq. (6). 

Appendix B: Random-Phase Approximation 

Here we consider the full spin and charge response 
functions taking into account the vertex correction due to 



various interactions. For d xz - and c? yz -orbital electrons, 
we consider the following on-site interaction 



V 



Hi = U E «*,rt n *,H + 17 E 



2 ^ 



S 7 > " Si 



J 1 



(Bl) 



where U and V correspond to the intra- and inter-orbital 
Coulomb repulsion respectively, J is the Hund's coupling, 
and J' is responsible for pair hopping. Note in (Bl), 
for each site i, rii r = rii^ + ni^ r ^. Within Random- 
Phase Approximation, irreducible spin response function 
in (Al) will be renormalized to be 

X.(q,«) = X?(q,w)[I ~ T s x° s (q,"T\ (B2) 

where I is the 2x2 unit matrix, x® is a 2 x 2 matrix 
formed by the intra- and inter-orbital spin susceptibilities 
defined in the last term in (Al), i.e., [Xs]rt = xirti an< ^ 
the vertex matrix 



r _ ( U J/2 

1 " ~ { J/2 U 



(B3) 



While the irreducible charge response function in (A2) 
will be renormalized to be 

X c (q, w) = x° (q, w) [I + T c x°(q, w)] " 2 , (B4) 

where [x%t = X°c, r t and 



U 2V 
2V U 



(B5) 



In the calculation of Imx s in Sec. Ill, we have chosen 
U = 1.8 and J = for simplicity. We do not actually 
calculate the RPA-corrected Im% c in the current context. 
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